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0> \ Abstract: The decays r -^ Svr + z/^ are analyzed at one loop order in the framework 

^ \ of Generalized Chiral Perturbation Theory, in order to test the sensitivity to the size 

I I of spontaneous chiral symmetry breaking parameters, contained in the S-wave. The 

Q^' latter, due to a kinematical suppression, at threshold, of the P-wave, is relatively 

i-^ ' large enough to be detectable at high energy machines, through azimuthal left-right 

asymmetries. This quantity (for the tt'-k^ti^ mode), integrated from threshold to 

/\' Q^ = 0.35 GeV^, varies from (17±3)% in the standard case of large condensate up to 

cd ' (40±5)% in the extreme case of tiny condensate. The feasibility of such measurement 

at high luminosity colliders {e.g. CLEO) is discussed. This method provides a 

completely independent cross-check of forthcoming experimental determination of 

the quark condensate, based on low energy vrvr scattering. 
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1. Introduction 

The mechanism of dynamical (chiral) symmetry breaking and the identification of 
its experimental signatures is one of the central subjects in modern particle theory. 
The problem is to understand non-perturbative alternatives to the spontaneous sym- 
metry breaking triggered by elementary weakly coupled scalar fields, whose nature is 
essentially perturbative. The problem arises at different energy scales characteristic 
of the Standard Model and its extensions. At low energy QCD, E < 1 GeV, the 
nonperturbative dynamical breaking of chiral symmetry (DBCHS) is unavoidable, 
since QCD does not contain elementary weakly coupled scalars. Yet, the mechanism 
of DBCHS, in particular the role and importance of the quark-antiquark condensate 
(qq), axe not understood theoretically and they remain poorely known experimen- 
tally [1]. At higher energy scales E ~ (10^ — 10^) GeV, a dynamical electroweak 
symmetry breaking is of interest as a possible alternative to an elementary weakly 
coupled light scalar Higgs protected by Supersymmetry This alternative assumes 
the existence of a new QCD-like (confining) gauge theory which breaks dynamically 
its chiral symmetry at a scale Agb ~ (2 — 3) TeV. The Goldstone bosons would then 
transmit this DBCHS to the electroweak sector by the Higgs mechanism (for a recent 
review see Ref. [2]) . In order to make such a scenario specific and to investigate its 



consistency with electroweak precision measurements at low energies {E <^ Agb), one 
usually takes the presumed mechanism of DBCHS in QCD as a prototype, assuming 
in addition, a simple dependence on the scale and on the number of (techni-) colors 
Nc and of light (techni- )fermion dublets A^d [3]- If the mechanism of DBCHS in 
Q CD-like theories turned out to be different from the general expectation (concern- 
ing e.g. the role of chiral condensates in a dynamical generation of fermion masses, 
the possible existence of new phase transitions as a function of A^c and A^d and their 
influence on the spectrum of bound states, . . .), the estimates of oblique and vertex 
corrections to the tree level Standard Model [3] could be substantially modified^. 
Below the scale Asb, at which new degrees of freedom ( "techni-hadrons" ) start to be 
visible, the precision tests of models of dynamical electroweak symmetry breaking 
can hardly be dissociated from the investigation of the mechanism of DBCHS in low 
energy QCD. The importance of dedicated experimental tests of the latter for the 
low-energy QCD phenomenology as it appears e.g. in the calculation of weak matrix 
elements, determination of CKM matrix elements or estimates of running light quark 
masses, has been often mentioned. (It has, for instance, been pointed out [5] that if 
the condensate (qq) does not dominate the expansion of symmetry breaking effects, 
the SM evaluation of e'/e could lead to a significantly higher value than in the case 
of the standard scenario of DBCHS). Actually, model independent tests of the role 
of (qq) in QCD are of more fundamental and larger interest: they may concern the 
mechanism of dynamical electroweak symmetry breaking at the TeV scale and the 
problem of the origin of masses. 

The purpose of this paper is to propose a new test of the importance of the quark 
condensate in the mechanism of chiral symmetry breaking in QCD. The experimental 
signature of this mechanism can be systematically investigated using the technique 
of the effective theory [6] , and its low energy expansion - xPT [7, 8] . The standard 
version of this expansion (S^PT) assumes a dominant role of {qq) that is a small 
deviation from the Gell-Mann-Oakes-Renner (GOR) relation [9]. If the condensate 
is less important the standard low energy expansion of the same effective Lagrangian 
has to be modified leading to the Generalized xPT (G^PT). The latter is the proper 
framework for measuring (qq) since it makes no assumptions on its size nor about 
the deviation from the GOR relation^. 

In this way the low energy tttt scattering in the S-wave has been identified as 
being particularly sensitive to (qq) [10, 11, 12]. The corresponding experimental 
tests now become possible, due to i) the control of the precision of xPT which now 
reaches two-loop accuracy and is able to produce a precise prediction of the standard 



^A remark in this sense has been recently made in Rcf. [4], concerning the value and the sign of 
the parameter Liq. 

^The quark condensate by itself is not an observable quantity. What can be actually measured 
is the renormalization group invariant product m{uu + dd)rh=o in units of F^M^, that is precisely 
the deviation from the GOR relation. 



scenario [7, 13] and on the other hand to interpret any substantial deviation from this 
prediction in terms of unexpectedely weak values of quark condensate [11]; ii) the 
fact that the needed precision (5 % or better for the S-wave scattering length) can 
be reached: in the on-going tt^ — n^ atoms lifetime measurement at CERN [14], in 
the new precise K^^ experiments, E865 at BNL [15] and, in the near future, at the 
Frascati factory Da0ne (with the KLOE detector) [16]. 

In view of the importance of this issue it is interesting to identify independent 
sources of similar experimental information. In this paper a proposal is elaborated 
in detail concerning the decays r ^ Svr + z/,- at low invariant mass of the hadronic 
system Q^ < 0.35 GeV^. The S-wave is proportional to the divergence of the axial 
current and it measures the importance of chiral symmetry breaking away from the 
pion pole. However, when considering this process in xPT, one is faced with two 
main difficulties: i) the smallness of the branching ratio in the threshold region, 
which is the only domain where xPT is applicable, and ii) the chiral suppression 
of the S-wave compared to the P-wave (which, moreover is enhanced by the ai reso- 
nance). The former problem now becomes less stringent, as pointed out in Ref. [17], 
thanks to the large statistics already accumulated in the present machines and to 
the improvements expected for the near future (for a review see Ref. [18]). As for the 
latter, it turns out (cfr. Ref. [19]) that the P-wave, although dynamically dominant, 
is kinematically suppressed, near the threshold, compared to the S-wave. These two 
facts together allow the S-wave contribution to show up in a clean and detectable 
way. The plan of the paper is as follows. In Section 2 the generalized SU(2)xSU(2) 
chiral lagrangian is constructed and renormalized at 0{p'^) level. In Section 3 we 
will recall the kinematics of the r -^ SnUr decays and the definitions of the structure 
functions, following Ref. [19]. In Section 4 we compute the hadronic matrix element 
at one-loop of G^PT, and in Section 5 we discuss the relation with tttt observables. 
The numerical results for the structure functions and their dependence on the quark 
condensate are shown in Section 6. In particular we will focus on the azimuthal 
asymmetries, which are very straightforward observables to extract experimentally. 
In Section 7 we discuss the relevance of our results for the QCD sum rules used to 
extract the quark masses. Finally we will summarize with some conlcuding remarks 
in Section 8. 

2. The generalized expansion of the SU(2)xSU(2) lagrangian 

XPT is an effective theory having as degrees of freedom only the lowest energy 
excitations of the spectrum of QCD, that is the Goldstone bosons. In the case 
of SU(2)xSU(2) chiral symmetry these are the pion fields tt^{x). They are collected 
in a unitary 2x2 matrix 

U{x)=exp "^J , J = 1,2,3, (2.1) 



with T^ the Pauh matrices, and F the pion decay constant in the chiral hmit. The 
external sources coupled to the vector, axial-vector, scalar and pseudoscalar quark 
currents are also represented by hermitian 2x2 matrices 
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^M = II^M^''' «M = II«M^^ s = J2s^t\ P = J2p't^- (2.2) 

i=o i=i j=o i=o 

The scalar source (which contains the quark mass matrix) and the pseudoscalar one 
are grouped in a complex matrix x = s + ip. The effective lagrangian is the most 
general one which is invariant under chiral symmetry and it contains an infinite 
number of terms; however, only a finite number of them will contribute at a given 
order of the low energy expansion. Since the latter is a simultaneous expansion in 
powers of momenta and quark masses, one has to establish a power counting rule 
for the quark masses. This can be done by considering the Taylor expansion of the 
pion mass squared in powers of quark mass (we will always stick to the isospin limit, 
rriu = iTi^ = m) 

Ml = 2Bm + AAm^ + ..., (2.3) 

which is valid up to chiral logarithms. The coefficient of the linear term is propor- 
tional to the quark condensate, 

5 = - lim ^ = - lim ^. 2.4 

»7i^0 F^ m^O F^ 

The standard assumption is that the latter is large enough for the first term to 
dominate the series (2.3), yielding the standard counting rule rh ~ 0{p'^). In general, 
however, nothing prevents B to be much smaller than one may wish; in this case the 
counting rule has to be modified and, to be consistent, one has to consider formally 
i? as a small parameter (notice that the parameter B has been removed from the 
definition of the matrix x, compared to the standard notation). The generalized 
counting is then defined by 

5 ~ m ~ X ~ 0{p), (2.5) 

so that the expansion of the effective lagrangian in G^PT reads 

r^ = £(2) + £(3) + £(^) + ... (2.6) 

with 

£W= J2 B^'^^Ki) (2-7) 

k+l+n=d 

and C(^k,i) containing k derivatives and / powers of quark masses. Defining the co- 
variant derivative as usual 

D^ = d^U-i [v^, U] - I {a^, U} , (2.8) 



and denoting by (. . .) the trace of the flavor 2x2 matrices, the leading, order 0{p^), 
lagrangian reads 

£(2) = ^ ^(D.U^D^U) + 2B{U^x + X^U) + A((f/tx)' + {x^uf) 

+ZP{U\ - x^UY + h, {xh) + K (det X + det x^) ] . (2.9) 

Terms which would be of order O(p^) in the standard counting are already present 
at the lowest order in the generalized case, reflecting the fact that the linear and 
the quadratic contribution in Eq. (2.3) can compete. Up to a given chiral order the 
SxPT lagrangian is a subset of the G^PT one. The last two terms are high energy 
counter terms, needed for the renormalization of the theory, but they will not appear 
in the expressions of matrix elements of pions on the mass shell. Using the equations 
of motion, field redefinitions, the identity 

M2-M(M) + detM = (2.10) 

valid for any 2x2 matrix M, and the fact that 

D^U^D^^U = ^{D^U^D^U), (2.11) 

to absorb the redundant terms, we arrive at the following expressions for the C^^^ 
lagrangian^ 

C^'^ = \F'{eHD,U^D^u{x^U + U^x))+P?\{x'uy + {u^xy) 

+P?\{x^u + u^x) x^x) + P?\x^u - u^x){{x^uy - {u^xf) 

+pf {x^U + U^x) ( (X^U) ' + (U^x) ') + P?^ ix^x) {X^U + U^x) } , (2.12) 
and for £(^) = £(4,0) + -^(2,2) + -^(0,4), 
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C^A,o) = jiD^U^D^Uf + j{D^U^D,U){D^U^D^U) + k{F;^,UF'^'>^^U^) 



- (2/^2 + ^/s) {Fj^.F'"'^'' + F/iz/^F^''^'^). (2.13) 

C-i2,2) = \f' {ar{D,U^D^U (xh + U^x^u)) + a2{D,U^Ux^ D^UU^x) 
+a,{D,U^U (x^D^x - D^xh) + D,UU^ (xD'^x^ - D^XX^)) 



■'The same lagrangian, up to and including 0{p^), has been found, independently, by Marc 
Knecht [20], to whom we are indebted for useful correspondence. Notice that in this reference 
different notations have been chosen for the pure source terms h'^ and ft-(2,2)- 



+h2{D^U^xD^U^X + X^D^Ux^D^U) 
+ ' 



, h{U^D^xU^D^x + D^x^UD^^x^U) 

+c^{D,U\ + x^D^U) {Df^U^x + X^D>^U) 
+C2{D,x^U + U^D,x){D^'U^x + X^D^'U) 
+cs{D,x^U + U^D,x){D^X^U + U^D^x) 
+c,{D,U^X - X^D,U){D^U^x - X^D^U) 
+c,{D,x^U - U^D,x){D^X^U - U^D^x) + h^2a){D,x^ D^^x)} ■ (2.14) 

£(0,4) = \F^[ei{{x^u)'+{u\)^) + e2{x\(x^Ux^U + U\U\)) 
+ez{xhU^XX^U) + /i((x^f/)' + (f/^x) V 

+f2{(x^u)' + (f/^x)')(x^f/ + f/^x) + h{xh (x^u + f/^x))(x^f/ + u\) 
+U{{x^uy + (f/^x)')(x^f/ + u^xf + M{x^uy - {u^xy){x^u - u^x) 

+hA{x^XX'^x) + hix'^x) (det x + det x"^) 

+/i6 (det X + det x"^) + /ly (det x - det x"^) }• (2.15) 



The lagrangian £(4,0) is the same as in Ref. [7]. The remaining terms of the 0(p*^) 
lagrangian of Ref. [7] , related to the exphcit chiral symmetry breaking sector of the 
theory, have been redefined; the correspondence reads: 

A = ^{k + h), Z^ = |J(-/3-/4 + /7), ^^'^ = ^h, (2.16) 

ho = -p:^{-k + /ii + /is), K = -p^ih + h- hs). (2.17) 

The divergences at 0(p'^) in dimensional regularization are absorbed by the following 
renormalization of the coupling constants*^ 

const = const*" + FconstA, (2.18) 

with 

and the Fconst? calculated with the standard heat kernel technique, are collected in 
the Table 1. We see from the table that, due to the chiral counting of B, the 
divergences come all at order p^, as they should. 



^A more detailed description of the generating functional, as well as the correspondence between 
the SU(2) and the SU(3) constants, will be given elsewhere [21]. 
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Table 1: The (3 function coefficients of the low energy constants. The scale dependence is 
fi (d/d^)const'" = -1/ (I67r2) Tconst- 

3. Kinematics 



There are two different charge modes in the decay r -^ SttUt-, the 2tt^tt mode and 
the all charged one, 27r^7r+. The relevant hadronic matrix elements are 



H^'-{pi,P2,P3) = {7i^{pi)n'{p2)7r-{ps)\A-^\0), 
H;'+{pi,p2,P3) = (7r'(pi)7r-(p2)7r+(p3)|^;|0), 



(3.1) 



with A = u'jfj^'-f^d. Both these matrix elements can be simply expressed as follows: 



H. +(pi,P2,P3) = V2\h+ \p3,P2.Pi) + H+ °(P3,Pi,P2) 



H'r{pi,P2,P3) 



H. '^iP2,P3,Pl)+H ^iPuP3,P2)- H +(pi,P2,P3) 



V2H;-\px,P2, 



P3) 



with 



and 



H^-%Pi,P2,P3) = {n^{pi)n-ip2)n%p,)\Al\0), 



^l = M7m75M - c?7m75C?- 



(3.2) 
(3.3) 

(3.4) 
(3.5) 



H^^^ contains directly as pole the invariant amplitude A{s\t,u) of the elastic tttt 
scattering. Both Eqs. (3.2)-(3.3) follow from isospin and charge conjugation sym- 
metry. The most general Lorentz structure of the matrix elements compatible with 
G-parity invariance is 

Hf{puP2.P,) = V,,F^'^{p,,P2,P3) + V2,Ff\p,,P2.p,) + K4M^r(Pl,P2,P3), (3.6) 

where hfs (hadronic final state) stands for 00— or h and 

vT = P^-p!^-^^g^ (3.7) 

v^ = pt+p'i+p''^ = Q^. 

The G-parity violating contribution ie^''^'^''^Pia.P2i3P?,-iF^{.PiiP2iPz) would be present 
if m^ 7^ TJifi. However, such a term could only originate from the WZW anomaly, 
which does not depend on quark masses. Therefore it will start to contribute only 
at the loop level, that is, in the case of the anomaly, at order 0{p^). Notice that, 
in general, m^ — m^ corrections are expected to be of smaller relative importance in 
the case of small condensate, since the ratio {m^ — m^) / {m^ + m^) should decrease. 
For both the 27i^7i^ and 27r°7r~ final states, Bose symmetry requires 

Ft{puP2,P3) = F^'%P2,PuP3). (3.8) 

Fi {F2) and F4 correspond respectively to a final hadronic state of total spin 1 and 0, 
as it is easily seen in the hadronic rest frame. The differential decay rate is given by 



(270^ 

2Mr 



dr(r ^ Sttz/,) = ^— ^ \M |M<I>4. (3.9) 



Ai is the matrix element of the electroweak interaction 

V2" "^ 



M = Kd^i^'^^f , (3.10) 



with 

L,^ = U^riPur, Si.j7M75Mr(Pr, S^), (3.11) 

and d$4 is the invariant phase space of four particles 

^^^ = 256 (2vr)i^ M,^ ^dad7dcos^d., d., . (3.12) 

dQh is the solid angle element of the hadronic system in the r rest frame and can be 
written, after integration over the azimuthal direction, as 

dnh = 27r(dcos^), (3.13) 



where 9 is then the angle between the laboratory and the hadronic system in the r 
rest frame, a, (3 and 7 are the Euler angles which describe the orientation of the 
hadronic system in the laboratory frame (for detailed definitions see Ref. [19], which 
we follow closely). The hadronic invariant variables are 

Sl = {P2 + Pzf, S2 = (pi + Psf, S3 = (pi + P2)^ (3.14) 

with Sl + S2 + S3 = Q^ + 3M^. 

As shown by Kuhn and Mirkes [19], the matrix element squared can be written 
in terms of 9 independent leptonic and hadronic real structure functions Lx and Wx 

\M\' = V^,^J:LxWx. (3.15) 

^ X 

All the angular dependence is contained in the functions L^'s while the Wx^s only 
depend on the hadronic invariant mass Q^ and on the Dalitz plot variables si and 
S2- Four of the hadronic structure functions correspond to the square of the spin 1 
part of the hadronic matrix element. 



Wa = {xl + a;2)|Fi|2 + {x^ + xl)\F2\^ + 2(xia;2 - xl)Re{FiF;), 
Wc = {xl - xi)|Fi|2 + {xl - xl)\F2\^ + 1{xxX2 + a;i)Re(FiF2*), 
Wb = 1 [xiXslFip - X2X3\F2\^ + X3{X2 - xi)Re(FiF2*)] , 
l^i, = -2a;3(xi+X2)Im(FiF2*), 



(3.16) 



one is the square of the spin component 

WsA = Q^\F4\\ (3.17) 

and the remaining ones are the interference between the spin and spin 1 components 

WsB = 2v^ [xiRe(FiF;) + X2Re{F2F:)] , 
Wsc = ~2./Q^[xilm{F,F:)+X2lm{F2F:)] , 



WsD = 2VQ^X3 [Re{F,F:) - Re(F2F;)] , ^^'^^^ 



WsE = -2VQ^X3 [Im{F,F:) - Im{F2F, 



4, 



The Xi are kinematical functions which are linear in the pion three-momenta in the 
hadronic rest frame^, xi = Vi, X2 = V2 and 3:3 = Vi: they vanish at the production 
threshold and can be expressed as follows: 

X,X2 -Xl = (QV2 - S3 - M^/2) + (S3 - si)(s3 - S2)/(4Q2), 

a^i + a;2 = -3/2 v^, 

0:3 = Vh, (3.19) 

-ho = AMI - (2M2 - Sl - S2f/Q\ 
h = [S1S2S3 - M'^{Q' - M,2)2] /(hoQ'). 



■''This system is oriented in such a way that the x-axis is along the direction of p^ and aU the 
pion fly in the x-y plane. 



We can see, from the explicit expressions for the Wx^s, that the purely spin 1 struc- 
ture functions are suppressed, at threshold, by two powers of the x^'s, the interfer- 
ences of spin 1 and spin by one power whereas the purely spin structure function 
is not suppressed at all. This fact, as already pointed out in the introduction, allows 
us to have, at threshold, a much larger sensitivity to the S-wave, and then to the size 
of chiral symmetry breaking measured by the divergence of the axial current. 

4. The hadronic matrix element in GxPT 

We have seen how the hadronic matrix elements for both the charge modes can be 
expressed in terms of H^~^{pi,p2,P3). It is then sufficient to compute this quan- 
tity: from the residue of the pion pole we will immediately read the tttt scattering 
amphtude. The G^PT result at one loop level is 
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2)m 
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327r2F,2 
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~327r2Ff 
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1 + i^^>m - 2(e(2)m)2 + (3a^ + Aa^ + Ah\ + 26^ + Ac\ - 2c2)m^ 
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FJ \2 
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isi - -M^ - 2Am'\ r(si) + -^ 
3 6^ y V 1; ^ 2887r2 



S2 - Si, 



(4.1) 



The scale dependence of the l.e.c.'s is absorbed by the corresponding dependence in 
the standard loop functions J^{s) and Af (s), defined e.g. in Ref. [8], and in the 
chiral logarithm 

1 / M^ \ 

A;.. = :r:— .llog^ + 1 , (4.2) 



327r2 



/i^ 



10 



so that the total result is independent of the xPT scale. A(s3|si,S2) is the (scale 
independent) off shell tttt scattering amplitude, with si + S2 + S3 = 3M^ + Q^, 

B / 4 \ M^ 
487r2Fi V 3/ V 3 



1 



487r2F4 



/.-^ 



si-2M2)' + (s2-2Mf ' 



M^ + 32Am2s3 - 24Am2M2 + 112^^771^1 J(s3 
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1 \^ 1 
Ml + 4Am2 - -si j + -(S3 - S2)(si 

M2+4Am=^-is2) +^(S3-Sl)(52 



4M2) 
4M!1 



J{s,] 

J{S2) \, (4.3) 



Putting Q2 = ]\^2 Qj^g recovers the tttt amplitude written in the same form as in 
Ref. [11] (dropping 0{p^) contributions). The expressions of the scale-independent 
parameters a, /3, li and I2 in terms of the low energy constants are given by the 
following equations: 

■K 

+ (lee-; + 4e^2 + 32/[ + 40/2^ + 8/3^ + 96/1) ^' 
^Aa^Mlrh' - ^^J^^ {?>Ml + 20 Am^) log ^, 



:Ml = 2Bm + 4A"m2 + {^p\ + p^ + 20pl + 2pl) m^ 



1 + 2^(2),r^ ^ (2a^ + a^ + 403 + 26^ - 26^) m^ 



(4.4) 



/i 



^7r2F2 



/3 = 1 + 2e'^'''m - 4^ 



{2),r, 



-(2)^ - 2 



4:Ml 



327r2F2 



10- 






2 (3a^ + 2a3 + Ah\ + 26^ + Ac\) m^ 



', Ml 

""'IF 



(4.5) 



(4.6) 



aMl^ = 2Bm + 16A'm^ - AMH^'^^'^m 

TV T^ 2 '' 

+ (81p^ +pl + lUpl + 2pl) m^ 

-8Ml {2bl - 2bl -as- Acl) m^ 

+ 16 (6Aa3 + 16e; + e^ + 32/[ + 34/2" + 2f; + 72 f I) m* 



Ml 



32n^F^ 

1 
'32n^F^ 



AMI + 204Am2 + 528 



A^m^\ Ml 



Mt + 88 Am^M: + 528 A^rh 



2-4 



(4.7) 
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(4.8) 
(4.9) 



The parameter a contains the main information about the size of the quark conden- 
sate: at tree level it varies between 1 and 4 if (qq) decreases from its standard value 
down to zero^, while f3 stays always close to 1. Solving Eq. (4.7) for A and Eq. (4.6) 
for ^^^', we can express the form factors in terms of a and p. The result for the 
27r°7r^ mode in terms of the form factors -Fi(pi,P2,P3) and -^4(^1,^2,^3) introduced 
in Section 3 is: 






+ ^[1(53 - 2Ml) + ^W^Ml + S2- 2si) + h,Q^ 



-AMI - 2Q' + -S2 + 4si 



3F2 

1 



s, + ^Ml{a-l] 



-^T 77.00- 



2Q2 



-s, + -M^{a + 2] 



J(S3) + 

^(si), 



1 



3Ff 



2,0,0; , ■ 

S2--M2(- + 4) 



^(^2) 



(4.10) 



iP + l] 



Q^-M^ 



(3 {S3- ^M^) + ^Ml 



+ 



al + 2a3 + 2hl + 4c; - 4c2 - 2 ' ^ ^^ ^^03 - c^) 



Q' 



m 



Ml 



96n^F^ 
1 
^967r2Fj 

4o 



I^q; — 1) lo£ 



/i2 



Mi 



S2 + S1- 2Ml 



a 



3 "" 6 g2 



-(2g"-3s2-3si + 6M^)- 



3(S2 + Si - 2M^) 
5 (S2 - si)' 



6 g2 



tQ"^ -S2-S1 + 2Ml , , ., 

^-2h- ^—^ ^(53 - 2Ml) 



Q' 



+^^g^ 



[S2 - Sif + (S2 + Si - 2MI){S2 + Si - 4M2) 

S2 + Si - 2M2 



+ TT7T^ < S2 - 4M2 + ^ i,,^ iM2(4 - a) 



12F2 

S2 - Si 



Q' 



Q' 



2s2-Ml{A + a) 



1 



Q^-Ml 



3(s2-^(a + 2)M2' 



(S3-Si)(s2-4M2 



^(S2) 



^Thc relation between {qq) and a at one-loop level will be revisited in another publication [22]. 
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si - S2 r 



s.-AMl + ''^\:^^" Ml{A-a) 



Q' 



Q' 



2si-Ml{A + a] 



Q'^-Ml 



3(si-^(a + 2)M2' 



1 



1 



2F2 \Q^ 
1 

F^ 



Q^-S2-Si + 2Ml 



S3-S2)(Si-4M2) 

4 



■Ksi) 



S3 + ^M2(a-l 



-M^(a - 1)S3 + MM -a^ - — a + — 



7 



32 



9 



9 



16 
"9 



■Ks,) 



g^^i^lis^^O-^)^^^-^^^)^ 



+ 



S2 - 2M^)' + (si - 2M^ 



2\2 



487r2F4 

where we have defined the scale independent quantities 

327r2 M2 

7i /^ 



(4.11) 



(4.12) 



and 7j are the corresponding 13 function coefficients collected in Table 1. It is easy 
to verify, taking the S^PT expressions for a and /?, 



a 



(3^ 



1 + 
1 + 



1 



327r2 
1 



', Ml 



1 



6/^' + 211 



Ml 



87r2 



log- 



Ml 

/i2 



211 



Ml 

TT 

If' 



(4.13) 
(4.14) 



that one recovers the results of Ref. [17]. Notice that the form factor F4 is propor- 
tional to the matrix element of the divergence of the axial current. One easily verifies 
that the expression (4.11) goes to zero in the chiral limit. These two facts, together 
with the global renormalization scale independence, provides non-trivial checks of 
our calculation. It is remarkable that in Eqs. (4.10) and (4.11) all the contributions 
from the C^^^ lagrangian are absorbed in the definitions of a, /3, F.,^ and M.,^. 

5. Connection with low-energy tttt observables 

The renormalization-group invariant parameters a and (3 are essentially the (on shell) 
TTTT amplitude A{s\t,u) and its slope at the symmetrical point s = t = u = A/SM^. 
For this reason one expects both a and j3 to be less sensitive to the higher order 
chiral corrections than, say, the S-wave scattering lengths. They represent the 0(p*^) 
GxPT truncation of observable quantities a and P defined, including 0{p^) accuracy. 
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in Ref. [11], 



A{s\t, u) = AKMSF{s\t, u; a, (3] Ai, A2, A3, A4) + O 



Ah/ 'Uh 



(5.1) 



where the function Akmsf is expicitly displayed. Both a and j3 will be measured 
in forthcoming high-precision low-energy vrvr experiments (Dirac at CERN, E865 at 
BNL, KLOE at Da<l>ne). The values of a and j3 are constrained by the behavior 
of the phase 5\ in the p region, using the Roy dispersion relations [23] (resulting 
in the so-called Morgan-Shaw universal curve [24]). The parameters Ai, . . . A4 have 
been determined in Ref. [25] from medium-energy experimental vrvr phases, using 
a set of rapidly converging sum rules derived from Roy dispersion relations. The 
results depend very little on the size of {qq). Comparing the expression (5.1) with 
the explicit two-loop calculation performed in S^PT [26] one finds the relationship 
between the parameters Aj's and the low energy constants of the Standard lagrangian 
[27], 

\ = \i{Wl2M,...), 1 = 1,.. A. (5.2) 

Notice that the parameters (5.2) do not depend at all on the critical standard 0{p'^) 
constant ^3, which controls the size of {qq), and the dependence on the constant l^ 
shows up only at 0{p^) order. For this reason Eq. (5.2) can be inverted for /i and 
I2, yielding the values [27] 

[^ = -0.37 ±1.96, [2 =4.17 ±0.47. (5.3) 

where the error bars come from the uncertainty on the parameters Aj's and the esti- 
mated uncertainty of higher chiral orders. Strictly speaking the values of Eq. (5.3) 
have been obtained within the framework of S^PT. However, for the reason men- 
tioned above, the >C(4,o) constants h and I2 should not be affected by the variation of 
the quark condensate. Improvements in the precision of the determination of /i and 
I2 are expected from forthcoming numerical solutions of Roy Equations [28]. 

Using the currently available data on 5^ — 5\ extracted from the last K'^^ exper- 
iment (Rosselet et al. [29]) the best fitted values for a and (3 are [11]: 

^exp_ 2.16 ±0.86, /5""P = 1.074 ±0.053, (5.4) 

where the error bars are merely experimental. They have to be compared with the 
S^PT 0{p^) expressions (4.13)-(4.14), which, using for l^, and /4 the central values 
quoted in Ref. [26] 

ll{Mp) = 0.82 X 10~^ lliMp) = 5.6 x 10"^ (5.5) 

give 

a'* = 1.06, /5"* = 1.095. (5.6) 
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Taking further into account the higher order contributions, estimated in Ref. [27], 
the standard values for a and /3 become^ 

a"* = 1.07 ±0.01, /3'^* = 1.105 ±0.015. (5.7) 

Our strategy to detect the sensitivity to the size of (gg), cross-checking the future 
determinations from tttt scattering, will be to plot the structure functions for three 
sets of values of a and /3, 

a^* = 1.07, /3^* = 1.105, 

^exp ^ 2.16, /?<=^P = 1.074, (5.8) 

a = 4, /5=1.16. 

They correspond respectively to the S^PT predictions, Eq. (5.7), the central values 
obtained from the best fit to the Rosselet et al. data [29], Eq. (5.4) and the extreme 
case of tiny condensate (a = 4 with the corresponding /9 = 1.16 inferred from the 
Morgan-Show universal curve). 

The other numerical input for the form factors (not entering the tttt amplitude) 
that we will use is the value of /e, from the two-loop analysis of the vector form factor 
of the pion [30], 

/6 = 16.0 ±0.9. (5.9) 

The contributions from the low energy constants of £(2,2) and i3(o,4) will be considered 
as sources of error, with zero central values and errors according to naive dimensional 
analysis estimates, that is 

1 1 , ^ 

C(2,2) = ± -y, C(o,4) = ± -^, (5.10) 

H H 

where Ah is the typical mass of the resonance exchanged, say Ah ~ 1 GeV. 



6. Structure functions and azimuthal asymmetries 

In this section we will give the numerical results for the structure functions integrated 
over the whole Dalitz plot, at fixed Q^, 

Wa,C,SA,SB,SC = / (^SidS2WA,C,SA,SB,Sc{Q'^,Sl,S2), 

wd,e,sd,se = I dsi ds2 sign(si - S2) Wd,e,sd,se{Q^, Si, S2), (6.1) 



"^Both the values (5.3) and the SxPT prediction (5.7) are compatible with the central values 
referred to as "Set 11" given in Ref. [13] without error bars. Notice that the errorbars given in 
Eq. (5.7) do not include any uncertainty due to I3 and I4. Eq. (5.5) is taken as our definition of the 
standard case. 
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where the sign ordering is needed due to the Bose symmetry. The hmits of integration 
of S2 for a given value of si are s^ < S2 < s^ where 



4 = 2Ml + 2E1E3 ± 2^ El - M^^Ei - Ml (6.2) 

and 

E, = ,^ , E, = -^ (6.3) 

are the energies of the particles 1 and 3 in the 2-3 rest frame {i.e. in the system 
where P2 + Ps = 0), while the bounds of si are 

AMl<.s^< ^Q^-Ml^ . (6.4) 

Due to the presence of two identical particles an additional factor 1/2 has to be 
understood in the integration over the Dalitz plot. The results for both the 

charge modes are shown in Figs. 1 and 2. In both cases wa and Wc are practically 
indistinguishable at the scale of the plots. The reason for that is, as already pointed 
out in Ref. [17], that 

Wa-Wc = 2xI\F^-F2\\ (6.5) 

where Fi — F2 is antisymmetric with respect to the exchange of the particles 1 and 2, 
and, as it is clear from Eq. (4.1), the antisymmetric part starts at order 0{p'^). Wa 
is the main contribution to the total decay rate: 



dr 



^F Kid 

128M^ {2-Ky 



Ml - Q' 



Ml + 2Q ,,,,^2^^^2 



3M2 



W{Q')dQ' (6.6) 



where 



3A/f2 



and we can see that most of the other structure functions are much smaller compared 
to it. The solid, dashed and dotted lines correspond to the three sets of values of 
Eq. (5.8) discussed above. We see that all the spin 1 structure functions practically 
do not depend on (gg) (except wu for the all charged mode, which, being zero at the 
tree level, is very small and very much affected by the theoretical uncertainties). On 
the contrary a rather strong dependence is observed not only in the purely spin 0, 
but also in the interference terms wsbi ■■■,u!se- In particular, among the latters, wsb 
is the largest one and its magnitude is comparable to wa- Clearly, wsb is the most 
appropriate structure function where to look for the S-wave contribution. From the 
experimental point of view, the sensitivity to the size of (qq) should therefore be 
expected in azimuthal angular asymmetries (see Ref. [32]). These are obtained by 
integrating the differential decay rate (3.9) over all the variables except Q^ and the 
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Figure 1: Integrated structure junctions j or the 211^ ii^ charge mode, wa andwc (here and 
in the next figure) are indistinguishable at the scale of the plot. Solid, dashed and dotted 
lines refer to the three sets of values for a and [3 of Eq. (5.8). They correspond respectively 
to the standard predictions, the central values extracted from the Kei experiment, and to 
the extreme case of tiny condensate. 

azimuthal angle 7, using the explicit angular dependence of the functions Lx found 
in Ref. [19]: 



dr 



n2 



128M^ (27r; 



5 ud 






^^^m^.)wm^Q^% (6. 



with W{Q'^) defined above and the azimuthal distribution, normalized to 1, 



/(Q2,7) = 1 + \2{Q^^) (C^LR cos 27 + C(jD sin 27) 
+\i{Q'^,(3r) (Clr cos 7 + Cud sin 7) . 



(6.9) 



The asymmetry coefficients C'-^^, C(jj-), Clr and Cud in Eq. (6.9) are related to the 
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Figure 2: Integrated structure functions for the 2it tt^ charge mode. The legend is as in 
the previous figure. 



Kuhn and Mirkes' structure functions by the relations' 
1 A Q^\ 3M2 wc 



^^^ 3 1^^ M^JM^ + 2Q^W 



Q2 \ 3M2 WD 



c =-Mi 



Ci 



LR 



n 3M^ wsB 
'4M2 + 2Q2 w '' 



a 



UD 



TV 3M^ WSD 

4 M2 + 2Q2 w ' 



(6.10) 



The interesting quantity is the left-right asymmetry, whose coefficient is related to 
WsB- The funtions Aj in Eq. (6.9) result from the integration over the r-decay angle, 



K{Q ) Pt 



1 dcos^ 
-1 2 



Pi (cos ^), 



where 13^ is the r velocity 



(3r 



1 



M2 



\ beam 



(6.11) 



(6.12) 



*We are neglecting the polarization of the t's, which is justified provided that the Z exchange 
can be neglected (far from the Z peak). 
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and the angle ip is defined in Ref. [19]. For r's produced at rest, which is the case 
for the r-charm factories, this angle is cosip = 1, so that the functions Aj are equal 
to 1, whereas in the ultrarelativistic case, jSr = 1, they become 



Ai(g',i) 



A2(g',i) 



1 



(M2 - g^ 



iM^-Q'' + 2M',QHog^ 



Q^ 

M? 



-2 + 3 



M'r + Q' . ,^2 



Ai(g^i). 



(6.13) 



M2 - g2 
These two functions are plotted for different values of (3r in Fig. 3. We see from the 
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Figure 3: The functions Ai and A2 for different values of the r velocity (3r- 

figure that, for high energy machines, the functions Aj decrease the sensisivity in the 
threshold region, but the effect is not so drastic, especially for Ai, which multiplies 
the coefficient Clr- In order to take full advantage of the whole statistics we can 
consider the integrated left-right asymmetry, 

iVR-iVL 



AUQ' 



iVR + iVL 



(6.14) 
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where N stands for the number of events from threshold up to a certain Q^ and the 
subscript L (R) refers to the events with 7 e [0,7r/2] U [3tt/2,2tt] (7 G [7r/2, 37r/2]). 
This asymmetry^ for both the charge mode is shown in Fig. 4. The sohd hues 
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Figure 4: The left right asymmetry, for the two charge modes, as function of Q^. The 
three curves, solid, dashed and dotted, inside the errorband correspond to the three different 
values of a and (5 in Eq. (5.8) (like in Figs. 1 and 2). The corresponding lines outside the 
hands are the results at tree level. 



correspond to the standard case, the dashed hues to the central "experimental" values 
of Eq. (5.4) and the dotted ones to the extreme case of tiny condensate [Eq. (5.8)]. 
The bands represent the theoretical uncertainties, coming from the errors on the low 
energy constants [Eqs. (5.3), (5.9) and (5.10)^°] and from a variation of the xPT 
renormalization scale /x, at which the estimate (5.10) is supposed to hold, between 
500 MeV and 1 GeV. Particular attention has been payed to the sensitivity of the 



^The results for the asymnietries refer to the case of ultrarelativistic r's, where the asymptotic 
expressions of Eq. (6.13) can be used. 
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In the estimation of the error from Eq. (5.10) we have taken m ~ 20 MeV. 
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structure functions to the two constants /i and /2- One sees from Table 2 that 
for all the asymmetries considered in this paper the uncertainty coming from these 
two constants is well under control. This would not be the case for the up-down 
asymmetries, related to the structure functions wd and ^5^1. The lines outside 



X 


A 


C 


D 


E 


SA 


SB 


SC 


SD 


SE 


1 dw""^- 
Wx dh 

1 a4- 

Wx dh 


0.02 
-0.02 


0.02 
-0.02 


0.01 
0.3 


0.02 
-0.04 


-0.03 
-0.06 


-0.004 
-0.05 


0.04 
0.009 


-0.02 
0.2 


-0.03 
-0.2 


1 dw^ ^ 
wx dh 
1 dw^ 

Wx dh 


0.02 
-0.02 


0.02 
-0.02 


8 
8 


0.005 
-0.02 


-0.03 
-0.07 


-0.007 
-0.05 


0.02 
-0.006 


0.6 
0.6 


-0.05 
-0.07 



Table 2: The dependence of the nine integrated structure functions wx on li and l2- The 
derivatives are computed at Q^ = 0.35GeV^. 

the bands in Fig. 4 correspond to the result at tree level. We see that the one-loop 
corrections to the asymmetry are rather important, not exceeding however 30 %. 
Assuming a geometrical behavior of this series, the two-loop effects are expected 
to remain in the range of 10 %, which at present stage should be considered as an 
additional systematic error to the one shown in Figs. 4 and 5. Notice however that 
a better control of this systematic error is conceivable: i) most of this uncertainty 
comes from the spin 1 part of the matrix element (see also Fig. 6), and consequently 
an inclusion of the resonance contribution which dominates in this channel might shed 
more light on the problem; ii) the strong final state interaction in the S-wave can be 
treated more precisely by the standard procedure of unitarization (cfr. Ref. [31]). 

In the case of r leptons produced close to threshold (r-charm factories), the above 
asymmetries are somewhat larger since in this case the kinematical functions Ai = 
A2 = 1 increase the contribution of the low-Q^ region to the integrated asymmetries. 
Taking as reference the value Q^ = (600 MeV)^, Alr for the all charged mode 
goes from (60 ± 6)%, for tiny condensate (a = 4), to (28 ±4)% for the standard 
case. We emphasize once more that the reason for such a big symmetry breaking 
effect of the S-wave is the kinematical suppression of the P-wave near the threshold. 
Another asymmetry where the same argument applies is A'lr, plotted in Fig. 5 (the 
legend is the same as in the previous figure), which corresponds to the coefficient 
Clp}^. ^Lf^ is defined analogously to Eq. (6.14) with L referring to the events with 
7 G [vr/4, 37r/4] U [57r/4, 77r/4] and R to the complementary interval. Eventhough this 
asymmetry is of less practical utility than Alr, it is remarkable that the separation 
of the three cases near the threshold is entirely due to the purely spin structure 
function, which is suppressed by two powers of the quark mass. Now that we 
have the predictions for the azimuthal asymmetries, it is useful to know how many 
events are expected in the near-threshold region. Using the experimental value of 
the r lifetime, xPT is able to predict the integrated branching ratio, that is the 
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Figure 5: The asymmetry ^lr i^r the two charge modes, as function of Q^. The legend 
is like in Fig. 4- 



number of events from threshold up to a certain Q^, divided by the total number of 
r produced. This quantity is plotted in Fig. 6. It turns out that it is very similar for 
both the charge modes. The uncertainties coming from the low energy constants and 
from a variation of the xPT renormalization scale between 500 MeV and 1 GeV are 
invisible at the scale of the plot. On the other hand, we see that the convergence of 
the xPT series is much worse than in the case of the asymmetries, as already pointed 
out (typically the tree level contribution is half of the total one-loop result; cfr. the 
discussion in Ref. [17]). We thus conclude that, in contrast to the asymmetries, the 
one-loop xPT predictions for the branching ratio are less reliable, and should not 
be used as a basis for a precision test of xPT. With this in mind, according to 
Fig. 6, the number of events with Q^ < (600 MeV)^ is expected to be of the order 
of 10~^ times the number of r's produced. On the experimental side, the largest 
statistics, up to now, has been collected by CLEO, with a total number of r pairs 
of ~ 10^, that is a hundred of events below (600 MeV)^ for each charge mode of the 
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Figure 6: The integrated branching ratio for the two charge modes. The solid, dashed 
and dotted lines correspond respectively to the three values of Eq. (5.8) for a and [3. The 
thicker lines are the one loop results and the thinner the corresponding tree level ones. 

decays r -^ SnUr- With such statistics, measuring Alr for the all charged mode, 
one should already be able to extract a, and then the quark condensate, with a 
precision comparable to the present one, summarized in Eq. (5.4). According to 
M. Perl (see Ref. [18]), a factor ~ 5 increase in statistics can reasonably be expected 
in the near future, thanks to the new high luminosity 5-factories. That could allow 
a much sharper measurement of a, thus providing an interesting and completely 
independent cross-check of future precise determinations from tttt scattering. 



7. The size of the spin spectral function 

The scalar structure function wsa{Q'^) is proportional to the three-pion component 
of the spin r spectral function, that is the spectral function of the divergence of 
the axial current. This quantity enters in the sum rule determinations of the light 
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quark mass rh. The two-point function relevant for such determinations is 

MQ') = tJd^xe''i-^0\T{d^ (ui.i^d) {x)d'' (rf7.75«) (0)} |0), (7.1) 

whose imaginary part is proportional to the spectral function 

-Im^5(Q') = p{Q') = ^ E(2^)'^'(Q -Pn) k^l^^ {dl.l^n) |0)f . (7.2) 

Using the Cauchy's theorem and the analytic properties of the two-point function, one 
can relate the integral of the spectral function along the positive x-axis, to a contour 
integral of ip5{Q'^) at high energy. Evaluating the latter with the help of perturbative 
QCD together with the OPE technique for the non perturbative contributions, one 
arrives to the class of sum rules 



JO 2iri J\Q^\=so 



|Q2|=So 

The three-pion component of the spin spectral function reads 



PsAQ') = ;t^ \<'AiQ')+wsI^iQ')] ■ (7.4) 



St fn2\ ' C^-^) 
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This is the hadronic phenomenological input that one has to insert, in the low energy 
part of the l.h.s. of Eq. (7.3) in order to evaluate the sum rule. Due to the smallness 
of the wsA structure functions, a direct experimental information on p3^ is lacking. 
Therefore a model for the three-pion continuum has to be used, which should include 
the observed resonances in this channel. Yet the overall normalization of the spectral 
function modelized in this way, remains unknown. It has been proposed in Refs. [33, 
34] to normalize the hadronic model to the low energy behavior predicted by standard 
xPT. However, as we have seen in the previous section, the low energy behavior of 
WsA depends very much on the size of the quark condensate. In Fig. 7 we plot the 
ratio 

D(Q,/3)/r)2\ _ Ps^ W 

^^ ^ " PtiQ' 

of the 3-pion component of the spectral function for a given value of a and f3 divided 
by the standard xPT prediction. The two bands correspond to the second and third 
set of values of Eq. (5.8); also shown are the tree level results. We see that the low 
energy spectral function increases by a factor of ~ 12 in extreme case of a = 4, 
and the output of the sum rule for rn^ can be expected to increase accordingly. In 
Ref. [32] the threshold behavior of the same ratio has been analyzed neglecting M^ 
in the integration over the Dalitz plot. The latter approximation is actually not very 
good near the threshold and overestimates the result: at tree level, in the extreme 
case of a = 4, the ratio R^°''^\9M^) was found to be 13.5. However, this effect is 
compensated by the increase due to the one-loop corrections. 
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Figure 7: The three-pion component of the spectral function psniQ'^) for a = 2.16, (3 = 
1.074 (solid line) and for a = 4, /? = 1.16 (dashed line), normalized to the standard xP'^ 
predictions. The lines outside the errorbands are the tree level results. 

8. Concluding remarks 

The understanding of the mechanism of DBCHS in QCD requires the knowledge 
of the size of the quark anti-quark condensate. G^PT is the appropriate tool to 
pin down this quantity from low energy experiments. The physical observable that 
will presumably allow to do that is the low energy vrvr scattering. In this paper the 
SU(2)xSU(2) generalized chiral lagrangian is constructed to 0(p^) and renormalized 
at one-loop level. We have presented an alternative way of extracting (qq) , from the 
r decays into three pions, performing an analysis at one-loop level of G^PT. The 
dependence on the condensate is contained in the S-wave of these decays, which is 
an effect of explicit chiral symmetry breaking. However, due to a kinematical sup- 
pression of the P-wave, compared to the S-wave, the latter shows up as a detectable 
effect near the threshold, through the azimuthal left-right asymmetries. At Q^ large 
enough to allow for an experimental determination (Q^ ~ 0.35GeV^) it turns out that 
the integrated left-right asymmetry is larger for the all charged mode (see Fig. 4). 
In the case of ultrarelativistic r's (relevant, e. g. , for CLEO) the predictions of 
G^PT for this asymmetry, integrated from threshold up to Q^ = 0.35GeV^, range 
from (17 ± 2) % (in the standard case) up to (40 ± 2) %, depending on the size 
of (qq). These errorbars do not take into account the uncertainty coming from the 
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higher chiral orders. We estimate the latter to be ~ 10% of the total result. There- 
fore a measurement significantly larger than 20% would already signal a departure 
from the standard picture of large condensate. Unfortunately the number of events 
expected in this region is rather small: the branching ratio integrated from threshold 
to Q^ = 0.35GeV^ is of the order of 10^^ for each charge mode. The presently avail- 
able statistics (10^ t^t^ pairs of CLEO) will allow to measure this asymmetry, but 
not to distinguish between the standard and the extreme generalized case. However 
improvements can be expected from both theoretical and experimental sides. The 
theoretical predictions can in principle be sharpened by an analysis "beyond one 
loop", which should include, in the spin 1 form factor, the effect of the resonance ai 
(cfr. Ref. [35]) and, in the spin form factor, the strong final state interactions of 
the S-wave (unitarization procedure [31]). From the experimental point of view an 
improvement in statistics by a factor 5^10 seem reachable. That would allow a real 
measurement af the parameter a, related to (qq), providing a completely independent 
cross-check of future precise determinations from tttt scattering. 
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